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Abstract 

In a study of the reaction e~ ^W~W^ with the DELPHI detector, the 
probabilities of the two W particles occurring in the joint polarisation 
states transverse-transverse (TT), longitudinal-transverse plus transverse- 
longitudinal {LT) and longitudinal-longitudinal {LL) have been determined 
using the final states WW-^lvqq {I = e,fi). The two-particle joint polari- 
sation probabilities, i.e. the spin density matrix elements ptt, Plt, Pll, are 
measured as functions of the W~ production angle, Oyy-, at an average reaction 
energy of 198.2 GeV. Averaged over all cos 9w-, the following joint probabilities 
are obtained: 



Ptt = (67 ±8)%, 
Plt = (30 ±8)%, 
PLL = (3 ±7)%. 

These results are in agreement with the Standard Model predictions of 63.0%, 
28.9% and 8.1%, respectively. The related polarisation cross- sections ctti ^lt 
and aLL are also presented. 
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1 Introduction 

In the study of the reaction e-e+^W-W+ at LEP2, the DELPHI Collaboration [1] and 
other LEP experiments [2,3] have established that on average ~ 25% of W particles are 
longitudinally polarised, as predicted by the Standard Model. The present study measures 
the joint two-particle spin density matrix elements which give the probabilities that both 
Ws are transversely polarised (WtWt), both are longitudinally polarised (WlWl) or that 
one W is transversely polarised while the associated W is longitudinal (WtWl+WlWt) ■ 
In what follows, these correlations will be referred to as TT, LL and LT, respectively. 
This is a more detailed test of the Standard Model prediction for the W polarisation 
than those previously pubhshed. Production of longitudinal l^s is of particular interest 
because they are associated with the breaking of the electroweak symmetry. This study 
tests the theoretical prediction of the correlations and in particular, that the correlation 
LL is suppressed relative to LT. 

The previously published measurements of the spin-dependent correlations between the 
W particles in the reaction e~e^^W~W^ are by the OPAL collaboration [4] and by the 
L3 collaboration [5]. OPAL evaluated the cross-sections cttt, ctlt and from their 
data at 189 GeV with a comparatively low statistics. Their results are in poor agreement 
with the Standard Model. L3 used the WW^luqq and WW^qq qq events from their 
full LEP2 data set to establish the correlation between the decay planes of the two Ws. 
The correlation was found to be consistent with the Standard Model prediction. 

The analysis presented in this paper uses only the events in which one W decays into 
an electron plus a neutrino or a muon plus neutrino, while the other W decays into two 
hadron jets. These "semi-leptonic" events are kinematically well constrained and they of- 
fer the best available data for any detailed analysis of this reaction. The r semi-leptonic 
events are excluded because the uncertainties in their identification cause problems in 
WW correlation measurements. The fully hadronic final state WW^qqqq is also ex- 
cluded because of the uncertainties in jet reconstruction: the charges of the hadron jets 
cannot be well measured and the particles from the four jets tend to overlap in the space of 
the detector, resulting in uncertainties in associations between the Ws and the measured 
jets. 

A complete description of the polarisation states of the produced W particles is given 
in terms of the two-particle joint spin density matrix PaiAi'A2A2') where Ai and A2 are 
the helicities of the W~ and W~^ respectively. In terms of the W production helicity 
amplitudes, -F^f^^' ^P^^ density matrix elements arc defined by 

2^lxX\X2 I A1A2 1 

The normalisation is such that the trace of the matrix is unity. The initial state helicity 
sum runs over = ±1/2 and the W particle helicities run over Aj, Aj = ±1, 0. 

The hehcities of W particles can be determined from their centre-of-mass decay distri- 
bution asymmetries. The above definition of the PaiAi'A2A2' elements can be put [6] into 
the following form which is model independent and is directly applicable to experimental 
data corrected for backgrounds and detection efficiences: 



2 



da 



(1) 



d{cos 9w-)d{cos 9l)d(f)ld{cos 9l)d(j)2 
da / 3 y ^ 




d{cos6w-) VStt/ 



Here, 9w- is the angle of the W~ production with respect to the e~ beam, 91 and ^2 
are the polar decay angles of the W~ and in their rest frames and 0*, 02 ^-^^ the 
corresponding azimuthal decay angles, as shown in figure 5. The functions Dx\' are the 
theoretical decay distributions of the W particles in the helicity states specified by the A 
indices. Precise definitions of the angles 9^ 2 and 4>\^2 ^^^1 of the functions Dx\i relevant 
to the present analysis are given in section 3. It should be noted that the cross-section 
formula (1) is model-independent regarding the WW production process. 

The single- particle spin density matrix pw', derived from the WW spin density matrix 
PaiAi'A2A2' by summation over one of the indices (1,2), gives information about the po- 
larisaton of one W regardless of the state of the other. All nine elements of the single- 
particle spin density matrix can be determined using the data from semi-leptonic events. 
This was done in [1] and [2,3], where only the electron and muon decays of one W were 
used as the analyser of the W polarisation. The hadronic decays of W particles were not 
used because of the reasons outlined above and because the analysing power of polari- 
sation in such decays is greatly reduced as the result of the severe practical difficulty to 
distinguish quark jets from anti-quark jets. 

In the present study, it would ideally be desirable to measure the complete 9x9 ma- 
trix PaiAi'A2A2'- As the result of the limitations in the polarisation information from the 
hadronic W decays, only a small part of the joint WW spin density matrix can be mea- 
sured. It is possible to measure 5 (out of 9) diagonal elements (PA1A1A2A2) pl^s 9 complex 
off-diagonal elements. Instead of this incomplete set of individual matrix elements, the 
following three linear combinations of matrix elements are considered in this paper: 



The quantities pll^ Ptt and plt are composed of the diagonal elements of the full matrix 
and they can be interpreted as probabilities of the joint, i.e. correlated, polarisation 
states of the two Ws. The elements Ptt and piT do not distinguish between the + and 
- transverse helicities, and also the polarisations of the states W^Wj^ and W^W^I are 
combined. This is a reduced set of information about the WW polarisations but it is 
nevertheless very useful. The elements ptti Plt and pn can be measured in semi-leptonic 
WW events because cancellations in the sums (2) imply that the incompleteness of the 
polarisation information in the hadronic W decays does not matter [6]. In section 3 it 
will be shown how p^r, Plt and pll can be measured directly from the data without 
recourse to the individual spin density matrix elements in (2). 



Pll — Poooo ■, 

Ptt = P++— + P—++ + P + P++++ , 

Plt — P++QO + P00++ + P — 00 + Poo — ■ 



(2) 



3 



2 The Experiment, Treatment of Data and Simula- 
tion 

2.1 The Experiment 

The DELPHI detector is described in detail in [7,8] and its configuration during the 
LEP2 runs is given in [9] . The reference [9] gives a complete description of the selection 
of WW events in DELPHI. The present analysis uses the data taken at centre-of-mass 
energies between 189 and 209 GeV. The data are grouped into three sets at average 
energies of 189 GeV, 200 GeV and 206 GeV. The total integrated luminosity is 520 pb~^, 
and the luminosity-weighted average energy of all data is 198.2 GeV. Jet reconstruction 
algorithms as well as electron and muon identification are also described in [9]. 

2.2 Selection of Data and Monte Carlo Simulation 

2.2.1 Data selection 

The initial selection procedure for the channels fiuqq and euqq is based on the typical 
topology of those events. As already mentioned, events from the ruqq channels are not 
included in this analysis and thus they are a part of the background. The starting values 
of the data cuts are those hsted in [9] : 

(i) Visible event energy > 40% of the nominal centre-of-mass energy; 

(ii) The event transverse energy > 45 GeV; 

(iii) The event must have at least one muon or one electron identified; 

(iv) Electron or muon candidate's momentum > 20 GeV/c; 

(v) Charged lepton track angle with respect to the beam direction > 20°; 

(vi) The total track multiphcity in each hadron jet > 3; 

(vii) Reconstructed W~ and W'^ masses > 50 GeV/c^} 

The precise values of these cuts, in particular those on the event transverse energy and 
the minimum particle multiplicity in jets, were varied slightly for data taken at the three 
different average e~e~^ energies. Three-constraint kinematic fits were then performed to 
the reaction e~e'^—>-W~W^ on the selected data samples, requiring both Ws to have the 
same mass consistent with 80.4 GeV/c^. Cuts on the probability distribution were 
then applied, with the value of the cut (in the region 0.5% -1%) determined from the 
distribution in each of the three data sets. The final sample, taken at all beam energies, 
consists of 800 electron plus 880 muon events. This sample is somewhat smaller than 
that reported in [9] because we require full functionality of all parts of DELPHI. 

Particle momenta and angles obtained from kinematic fitting have been used in the 
analysis of this experiment. 

2.2.2 Event Simulation 

Simulation of events plays a crucial role in the experimental procedure to separate events 
corresponding to the reaction e~e+— the "signal", from backgrounds. The signal 
is defined by the three CC03 Feynman diagrams shown in figure 1. These account only 
for a part of the four-fermion processes contributing to the data. In this experiment a 
customised version [10] of the WPHACT [11,12] generator program was used to simulate 

^The masses are calculated from the results of a preliminary 1-constraint kinematic fit to the reaction e~e+— W+. 
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all the four-fermion processes. The DELPHI WPHACT program includes reweighting 
for the Double Pole Approximation (DPA) radiative corrections and the possibility to 
compute the matrix elements of different subsets of Feynman diagrams. The weights 
are the ratios of the squared matrix element for WW production only via the CC03 
diagrams to that for production via the full set of four-fermion processes. It is thus 
possible to simulate CC03 events corresponding to production via the CC03 diagrams 
with or without inclusion of other four-fermion processes. 

In addition to the four-fermion background, there is a significant two-fermion background, 
mostly from qq^y final states. This background has been simulated using the KK2F 
generator [13]. 

The generators were interfaced to the PYTHIA [14-16] hadronisation program. Large 
simulated samples, of the order of 10^ events, were produced by the programs listed 
above, interfaced to the DELPHI detector simulation program DELSIM [7,8]. 



Figure 1: CC03 diagrams, i.e. the lowest order contributions to the amplitude for W W^ 
production. (fi,2,3,4 stand for the appropriate fermions.) 



2.3 Reconstruction of events 

Event reconstruction efficiencies were determined using events simulated with the 
WPHACT Monte Carlo program adapted for DELPHI [10]. The efficiency is defined 
as the ratio of the number of reconstructed CC03 WPHACT Monte Carlo events to the 
number of generated events, including all four-fermion channels. 

Efficiencies are determined as functions of cos^vv-, cos^J and cos defined in section 1. 
The resulting reconstruction efficiency table is a 5 x 5 x 7 matrix corresponding to 7 bins of 

cos^viz-i as defined in section 4.2, and 5 equal size bins for each of cos9\ and cos The 
reconstruction efficiency matrices were determined separately for the electron and muon 
event samples at the three centre-of-mass energies. Typical values of the reconstruction 
efficiency, e, as a function of cos^^f- ^'^^ 0.65 to 0.74 when integrated over all cos^* and 
cos ^2- The values of e as a function of cos6\ (or of cos ^2)) averaged over cos^^-, vary 
over the range 0.6 to 0.8. 





e 



e 
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2.4 Treatment of Backgrounds 

The dominant backgrounds to the selection of eujj and jJ-i'jj events (where j imphes a 
hadronic jet) can be divided into two groups: 

(a) Events which result from problems in reconstruction or selection procedures. The 
dominant contribution to this class of events comes from two-fermion final states, in 
particular from qq'y. Other contributions come from neutral current four-fermion final 
states which might be misidentified as ei/jj or, more rarely, as iii^jj. Misidentified charged 
current events Ti/jj are also potentially a background in this experiment. 

(b) The non-CC03 charged current four-fermion contributions to the global fxvjj or 
evjj final states. In contrast to the background (a), this background would exist even 
in a perfect detector and an ideal experimental event selection procedure. A number of 
the charged current four-fermion diagrams in this category can interfere with the CC03 
amplitude^. 

Following the treatment in [9], the backgrounds requiring special attention are: 

(i) the r events, WW^rujj (belonging to group (a)); 

(ii) various four-fermion processes which for experimental or other reasons lead to the 
same final state as the reaction e~e'^^W~W'^ (such backgrounds may arise from groups 
(a) and (b)); 

(iii) two-fermion events, mostly gg7 interactions (group (a)). 
Discussions of each of these background sources follows in turn. 

• The background from source (i) has been investigated by passing the simulated Tvjj 
events through the normal analysis chain and requiring evqq or ixvqq fits. The background 
from this source turns out to be negligible as the result of kinematic cuts and fitting. 

• The background (ii) is due both to charged current and neutral current events. 

The charged current processes consist of the three CC03 diagrams plus seven diagrams 
with s-channel exchange of ^°/7, leading to the production of only one W . There are 
also ten charged current diagrams corresponding to t-channel processes which give rise to 
one W. The latter diagrams can give rise to backgrounds only in the electron final state 

The neutral current four-fermion states have two quarks and two leptons of the same 
flavour. If one of the electrons is not identified in the detector, the event may be classified 
as belonging to the channel WW — > evjj and may satisfy the criteria for an acceptable 
kinematic fit. Altogether, the four-fermion backgrounds affect the electron channel, euqq, 
more than the muon channel ixvqq. 

It was found that the non-CC03 four-fermion background in the real events could be 
efficiently reduced to a level less than 4% by tuning the kinematic cuts and the cuts, as 
described in section 2.2. The effectiveness of the removal of this class of background events 
can be demonstrated in the following way: Starting from a large sample of generated 
events, two data samples were produced making use of the event reweighting facilities in 
WPHACT. Sample A contained predominantly the CC03 events and sample B the non- 
CC03 four-fermion background. Each sample was processed through the experimental 
procedure described in section 2.2. The event ratio II = B/A represents the proportion 

*^As a consequence of the finite W width, any Monte Carlo generator of e~ e^—*W~W'^ events at low energy must be 
a four-fermion generator in order to satisfy gauge invariance. 
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of the four-fermion background in the WW signal to be expected in the final sample of 
the real data. This quantity is of the order of 3% and is weakly dependent on cos 9w- ■ 
A plot of n for WPHACT data at 200 GeV is shown in figure 2(a). The results for other 
run energies are similar. 

A further test of the effectiveness of the method for deahng with the non-CC03 four- 
fermion background is to apply the analysis to a simulated data set where the expected 
result is known. For this purpose, the element poo of the single-H^ spin density matrix p\y 
is evaluated. This test was carried out for data generated at all three run energies, but here 
only the results from WPHACT at 200 GeV are shown. In figure 2(b), triangle symbols 
are used to plot the value of poo obtained from the WPHACT generated events using 
all the four-fermion diagrams, with no cuts except for the final kinematic selection. 
Star symbols are used to plot the value of poo after the same generated events have 
been passed through the complete selection procedure described in section 2.2. The 
smooth sohd curve is from an analytic calculation using only the CC03 diagrams. The 
conclusion drawn from this is that the procedure adopted for analysing the data removes 
essentially all the four-fermion background, leaving events which are attributable to the 
CC03 signal.^ 

The residue of the non-VrVr four-fermion background in the real data at all run energies 
is estimated to be at the level of 3 ±2%, the uncertainty is due to combining results from 
different energies. 

• The background process (iii) is potentially very serious because the cross-section for 
qq'j production is about an order of magnitude larger than that for the WW signal. Al- 
though the topology of the qqj events is quite different from that of the WW events, 
the reconstructed events of the background can resemble and fit the 1^1^ reaction. The 
problem of how to suppress this background has been investigated using the KK2F Monte 
Carlo generator [13]. It is found that this background shows a very characteristic kine- 
matic signature in the distribution of the quantities cosG^ and cosGj^, where 0^ is the 
angle between the momentum vectors of the two leptons (one charged and a neutrino) 
coming from one W and Qjj is the angle between the jets from the accompanying W, 
with quantities defined in the laboratory frame (i.e. the reaction centre-of-mass).^ The 
presence of the qqj background shows up as accumulations of events at cos0// ~ — 1 and 
cosQ jj ~ — 1. 

After the usual cuts and kinematic fitting, the real data samples show peaks in the cos 0^ 
and cos Qjj distributions indicating small but non- negligible contamination from the qq'j 
background. Below it will be shown that the contamination is of the order of 10% at 
189 GeV. This contamination has to be evaluated accurately because it affects the angular 
distribution of W decays. The kinematics of the W decay in the laboratory frame of the 
reaction is such that the distributions of the quantities cosQu and cosQjj have large 
discriminating power against the qq^ background. This enables the magnitude of the 
background to be determined and also provides a means of reducing the background by 
applying cuts on these distributions or by applying suitable event weights. Only the 
method of weights has been used in dealing with this background. 

^The plot made with triangle symbols in Fig. 2(b) shows a strong dip near cos 9^- = — 1 and some fluctuations up to 

about cosOin,- = . These effects are caused by some of the four-fermion backgrounds in the simulated raw data. In 
particular, the background events whose W particles arc not genuine spin 1 states tend to interfere destructively with the 
true Ws in the poo evaluation. 

"^The momentum vectors are taken from the constrained kinematic fit. 
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Figure 2: (a) Plot of the quantity 11, defined in the text, as a function of 0086*1^- for data 
simulated at 200 GeV. 

(b) Plots of the density matrix element poo evaluated from four-fermion events 
generated at 200 GeV with WPHACT before (triangle symbols) and after (star symbols) 
the event reconstruction procedure, including application of data cuts and kinematic 
fitting. The smooth curve is from an analytic calculation using CC03 diagrams. 



Figures 3(a) and 3(b) show, respectively, simulations of the cosG^^ and cos9jj distribu- 
tions at 189 GeV. The curves labelled WW were obtained from the WPHACT simulation 
of CC03, while the curves labelled gg7 are from the KK2F simulation of this background. 
At this stage of the analysis, the relative magnitudes of the WW and gg7 components 
are still unknown; therefore, the two curves are normalised to the same number of events. 
Figures 4(a) and 4(b) show the results of a least squares fit to the relative contributions 
from WW and gg7, using the distributions from the simulated events, to the real data 
distributions. The fit to the 189 GeV data requires contributions of 90% from WW and 
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Figure 3: Distributions (a) of cosG^; and (b) of cos9jj in the gg7 background and in 
the reaction e~e^^W~W~^ from data simulated at 189 GeV. (The angles Qu and Qjj 
are defined in the text.) The curves labelled qq'j and WW are normalised to the same 
number of events. 

10 ± 2% from qq'-f, while the same analysis applied to the data at 200 GeV and 206 GeV 
require a 5 ± 2% qq'y background contribution at both energies.. 

This determination of the relative magnitude of the qq'-f background with respect to the 
WW signal enables a combined WPHACT plus KK2F simulation of the real data. The 
combined simulation was subsequently used to derive weighting factors, p, defining the 
purity of the WW signal in each bin of cos 6w-, cos 61 and cos ^2 5 using the same binning 
as that used for the reconstruction efficiencies, e. The application of the weight factors p 
is described in the following section. 
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Figure 4: Distributions (a) of cos 0^^ and (b) of cos Qjj for the real data at 189 GeV. 
The superimposed curves are the resuh of a fit using simulated samples of qq'y and WW 
events. 



3 Analysis 

W decays are well described by the V-A theory of the charged current weak interactions. 
The theory gives the functions D\x/{6*, 0*) required for evaluating the spin density matrix 
elements PaiAi'A2A2' by applying equation (1) to the data. The W~ production angle, 9w-, 
and the W~'~^ decay angles ^* 2, 0* 25 which specify the direction of the final state fermion 
in the rest frame of the W~ and of the final state anti-fermion in the rest frame of the 
W~^, are defined in figure 5. 

As already pointed out, fermions can be distinguished from anti-fermions in the W lep- 
tonic decays but not so in the hadronic decays. However, some of the decay functions 
Dxx'{0*,4>*) ell's invariant under the transformation which rotates the vector momentum 
of a fermion in the W rest frame into the direction of its opposite anti-fermion vector. 
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Figure 5: Diagram of the momentum vectors in the reaction plane of the process 
e~e~^^W~W~^ . The two right-handed sets of orthogonal axes {xi,yi,zi) and {x2,y2,Z2) in 
the respective rest frames of the W~ and are also shown. The polar angles 61 and 
62 of the vector momenta of the two-body decays of the W~ and are measured with 
respect to the axes Zi and Z2- 



These functions will be called symmetric and designated by D^'^\ Consequently, the sym- 
metric D functions and the symmetric parts of the non-symmetric D functions can be 
used to analyse the polarisation of the Ws decaying into the purely hadronic final states. 
The polarisation information obtained is thus incomplete but, nevertheless, it is useful 
and, in particular, the quantities ptt, Plt and pll can be obtained from the data using 
only the symmetric decay distributions in both associated Ws. 

The theoretical formalism for extracting ptt, Plt and pll from the data is based on 
equation (1) with two modifications: 

(i) Equation (1) is integrated over the full range of (pl and 02- This removes the 
functions Dxy having A' 7^ A and eliminates all non-diagonal elements of the matrix 
PAiA'jA2Ai,- The following three decay functions remain: 

= ^(i-cosr)2, 

Dooie*) = sin^r, (3) 

= ^(i + cosr)2. 

(ii) Furthermore, only the symmetric parts, -D^A' '^^ ^^e functions in (3) are to be used: 
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Pursuing the above steps, the following modified form [6] of equation (1) is obtained: 



+PLT{cosew-)\D\^^{ei)D^;\ 



^3^ 
v4. 



y(i(cos 



(i(cos ^i4^-)c?(cos 6'|)(i(cos 6'2) 



{PTT{co.ew-)Df^{ei)D 



(4) 



where ptt, Pll and p^r have been defined in (2). The above equation assumes CP in- 
variance, i.e. that the joint polarisations WJ^W^ and W^W^ are equal. The superscript 
(s, s) indicates that only the symmetric decays of both Ws are considered. 
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Figure 6: Angular distribution of W~ production. Points with error bars are the ex- 
perimental data; the errors shown are statistical only. The data binning is explained in 
section 4.2. The small circles show the predictions of the WPHACT simulation. 



The angular distribution of the W~ production is not used explicitly in the data analysis 
for measuring the correlations ptt, Plt and pll- However, it is appropriate to examine 
the shape of the distribution dN / d{cos6w-) as a test of the quality of the data. Figure 6 
shows the angular distribution of the W~ in the present data at all energies combined, 
compared to the WPHACT prediction. The agreement with WPHACT is satisfactory. 
The distribution shown also agrees well with that in the DELPHI publication on the 
WW production cross-section [9]. 

Before using the functions (3) to determine the intensities of different helicity states in the 
data, these functions must be transformed to a related set of functions with the property 
that each one is non-orthogonal to only one in the group of distribution functions (3) 
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and is orthogonal to the other two. These functions are called projectors [17] and they 
can easily be worked out from (3). The projector functions needed to extract p^r, Plt 
and pll are: 

Kl = 2-5cos^9*, (5) 
Ay = 5cos^r-l. 

These projectors arc normalised to give the spin density matrices in the standard represen- 
tation [6,17]. The quantities pll, ptt and p^r are obtained from the data by evaluating 
the following sums: 



-^EMoiMMo;,), (6) 



Pll = 
Ptt — 
Plt = 

where summations are over all events i. 
and the event weights w are defined as 



— Y,i Wi is the sum of all event weights Wj, 



P 

w = -, 
e 

where p is the WW purity factor, defined as in section 2.4, for each bin of 
{cos 9w-) cos 91, cos ^2)) ^-^d ^ is the event reconstruction efficiency in that bin, defined in 
section 2.3. In these formulae, index 1 refers to W~ and index 2 refers to W^. 

Monte Carlo studies have shown that the correlations p^r, Pll, Ptt extracted from small 
samples of data (such as we have at each of the three energy points considered) are 

subject to large statistical fluctuations. These fluctuations are much larger than those 
encountered in the determination of the single-W spin density matrix elements. Because 
of that, all 1680 semi-leptonic electron and muon events have been taken as one sample 
for measuring p^^, p^^ and ptt- The sum of weights of these events is 2844. 

The fact that events which we analyse here as one sample come from a spread of centre- 
of-mass energies presents no difficulty because the theoretical predictions which we want 
to test can be modified to take into account the spread. In particular, the Standard 
Model [6] predicts a negligible variation of ptt, Plt and pll over the energy range of this 
experiment. 



4 Results 

4.1 Systematic Effects and Errors 

(a) Data cuts. 

Systematic effects resulting from residual backgrounds in real events have been estimated 
by processing the data several times with small variations in the cuts and, separately, with 
various cuts. Variations in the resulting values of the spin density matrix elements plt-i 
Pll, Ptt amount to 10% of their statistical error. This is interpreted as the magnitude 
of the systematic uncertainty and it is neglected. 
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(b) Hadron jet reconstruction. 

Problems in hadron jet reconstruction [9] can give rise to a shift in cosOw- with the 
further consequence that the reconstruction efficiency is read from the wrong cell of the e 
matrix. This migration and its effects have been examined using simulated events. The 
resulting uncertainty on 9y\/~ is small: it varies from ±2° at large angles to ±5° at small 
angles, i.e. in the forward direction^. This is negligible by comparison with the sizes of 
the cos 9w- bins. The effect of this migration on the spin density matrix elements has 
been examined by moving the simulated events randomly by the above uncertainty in 
cosOiY-. The effect of this variation on the joint spin density matrix elements is 5% - 
8% of the statistical uncertainty and is therefore considered to be negligible. 

Also, jet reconstruction problems can produce wrong momentum vectors of the W 
hadronic decay products. This has been investigated by processing the same events 
using different jet algorithms as described in reference [9]. No statistically significant 
effect was found when comparing the spin density matrix elements distributions obtained 
in these tests. 

(c) Lepton charge determination. 

Tracks at small angle with respect to the e~ and e"*" beams are susceptible to wrong 
charge determination. This problem is essentially eliminated [1] by the 20° cut on the 
lepton polar angle (see section 2.2). The effect of this cut on the spin density matrix 
elements has been examined by simulation of events and was found to be negligible. 

(d) Radiative corrections. 

The effect of the initial state radiation is essentially removed by the appropriate en- 
ergy cut. The final state radiative corrections are implemented through the Double Pole 
Approximation in the WPHACT reweighting. The uncertainty due to the radiative cor- 
rections on the spin density matrix elements is negligible. This has been established 
by a comparison of the spin density matrix elements evaluated from events generated 
with WPHACT, including the corrections, with the same elements calculated analyti- 
cally [6,17] without the radiative corrections. This can be understood because the spin 
density matrix elements are ratios of quantities which are similarly affected by the radia- 
tive effects. 

(e) Use of a fixed W mass in kinematic fitting of events. 

Three-constraint kinematic fits of the reaction events are needed in order to separate the 
signal from a large background. Fixing the masses of both Ws in the reaction to the 
same value, as stated in section 2.2.1, does not cause noticeable distortions of angular 
distributions and other quantities needed for the physics analysis. This has been checked 
by comparing the results obtained from the three-constraint fits with those from the one- 
constraint fits of the same events. As a further test, somewhat different fixed values of 
the W mass were tried. Statistically insignificant differences in the results were found. 

4.2 Presentation of Results 

The measured values of ptt, Plt and pll as functions of cos^^- are shown in figure 7. 
Because the number of events in the negative hemisphere of W~ production is much 
smaller than that in the positive hemisphere, the data have been divided into two bins 
in the negative hemisphere and five bins in the positive hemisphere. (The positive hemi- 

^This uncertainty includes a problem in the simulation of charged particle tracks in the forward region of DELPHI. This 
problem and its solution are discussed in two recent DELPHI papers [18,19]. 
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sphere is in the direction of the e beam.) It is easy to derive from formulas (5) and (6) 
that the condition 

Pll + Ptt + Plt = 1 (7) 

is vahd on an event by event basis and is hence automatically satisfied by all data samples. 

The curves shown in the plots are the Standard Model calculations based on the CC03 
diagrams evaluated at 198.2 GeV using the expressions in [17]. The error bars shown 
in the plots are statistical, the systematic errors being negligible by comparison. Errors 
on all measured quantities are evaluated from the data as standard deviations. The 
distributions of the errors are approximately Gaussian. 

Some events can contribute negative numbers to the sums shown in (6) while the final 
result is positive. However, if a particular correlation quantity p is very small, the mea- 
surement errors, which have a Gaussian distribution, can lead to an overall result which 
is negative. This happens in three out of the seven measured values of pll presented 
here. They are very small negative quantities, consistent with zero within the measure- 
ment errors, i.e. | pll |< ^Pll, where 5pll is the measurement error. These negative 
values ol Pll are included in figure 7. The condition (7) still holds in these cases as a 
result of the properties of the projector functions (5). Of course, the physical quantities 
p should satisfy the condition p > 0. Several methods of adjusting the results to satisfy 
the positivity condition were tried, but not used in the end because they introduce biases 
that are more harmful than the small negative numbers among the results. 

The data bins in figure 7 are too wide to show the possibly rapid variations of the measured 
quantities very near to cos 9w- = ±1- However, the exact values of these correlations 
at cosOw. ~ ±1 follow from the conservation of angular momentum. Neglecting the 
electron mass, the vector and axial-vector interactions involved at the e~ and vertices 
cause the initial e~ and to interact only when their helicities are opposite. Thus, the 
initial system has total helicity ±1 and when the final state is collinear with the e~e^ 
beams its total helicity must be the same. This means that at d-^- = or tt, plt must 
be 1 and, at the same time, Pll — Ptt — 0. These values have been indicated with 
diamond-shaped symbols in figure 7. 

Because of the low statistics in this experiment, it is useful to examine the values of p^r, 
pLT and Pll averaged over all bins of cosOw- and to compare them with the predictions 
of the Standard Model. The results are shown in Table 1. 



p 


Measured fraction 


Standard Model 


Ptt 


67 ± 8% 


63.0% 


Plt 


30 ± 8% 


28.9% 


Pll 


3 ± 7% 


8.1% 



Table 1: Measured values of ptt-, Plt and pll-, averaged over cos6'iy-, compared with 
the predictions of the Standard Model. Errors are statistical. 



The differential polarisation cross-sections daTT/d{cos9\Y-), daLT/d{cos9w-) and 
daLL/d{cos9w-) are related to ptt, Plt, Pll through the equation 
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Figure 7: Variation of the spin density matrix elements TT, LT and LL with the cosine 
of the W~ production angle, Ow-^ at the average energy of 198.2 GeV. Points with error 
bars are the measured data, and the curves are the Standard Model CC03 predictions. 
The diamond symbols indicate the points predicted by angular momentum conservation, 
as explained in the text in section 4.2. 



~T( a — ^ = ^7? a — rPTTicosew-), (8) 

plus the analogous expressions for the states LT and LL. The first term on the right-hand 
side of (8) can be replaced by the data points from the angular distribution shown in 
figure 6, normalised to the total cross-section for e~e^^W~W^ at the average energy 
of the experiment. The relevant measured cross-section is a = (17.07 ± 0.57) pb at 
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Figure 8: Variation of the differential cross-sections TT, LT and LL as functions of the 
cosine of the W~ production angle, Oiy- . Points with error bars are the measured data; 
the curves are the Standard Model CC03 predictions for the average energy of 198.2 GeV. 



198.2 GeV. It has been obtained by an interpolation of the DELPHI measurements [9] 
which cover the range 161 to 209 GeV. With this result and the measured ptt, Plt, Pll 
as functions of cos6]y-, one obtains the differential polarisation cross-sections shown in 
figure 8. Integration over cosOw- yields the total polarisation cross- sect ions (Jtti ^^lt 
and shown in Table 2. 

These results are in good agreement with the Standard Model predictions. The polar- 
isation fractions expressed in terms of the density matrix elements (Table 1) and those 
expressed in terms of the cross-section ratios (Table 2) are different expressions of the 
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two-particle polarisation correlations in the reaction. Measurement errors of the polarisa- 
tions are themselves strongly correlated. The extent of all the correlations will be shown 
in figure 9(b). 



a 


Measured Cross-Section 


Standard Model 


(JTT 


(12 ± 1) pb 


(10.57 ± 0.05) pb 




( 4 ± 1) pb 


(4.95 ± 0.02) pb 




( 1 ± 1) pb 


(1.40 ± 0.01) pb 


aRatios 


Measured Ratio 


Standard Model 


<7tt/ C" 


0.70 ± 0.06 


0.625 ± 0.003 




0.23 ± 0.06 


0.292 ± 0.001 




0.06 ± 0.06 


0.083 ± 0.001 



Table 2: Measured values of the total cross-sections cttt, ctlt and an, at the average 
energy of 198.2 GeV, compared with the predictions of the Standard Model. Also: Mea- 
sured values of the ratios cttt/c", c"lt/c" and aiL/a, at the average energy of 198.2 GeV, 
compared with the predictions of the Standard Model. 



5 Discussion and Conclusions 

The measurements of the total cross-section for the reaction e'e'^^W^W'^ by ALEPH, 
DELPHI, L3 and OPAL [20] , have successfully tested the predicted gauge coupling cancel- 
lations in this reaction. The direct measurement of the longitudinal cross-section for the 
single (i.e. uncorrelated) W production [1-3] is an important test of the same cancellation 
because the cause of the potential divergence of the cross-section is in the longitudinal 
parts. The measurement of the joint WW hehcity states presented in this paper is an 
advance towards a more complete test of the Standard Model in the context of the WW 
reaction. Our measurement of a^e^e^^WLWi) is consistent both with zero and with 
the small value predicted by the Standard Model. This confirms that the probability of 
longitudinal W production is predominantly in association with a transverse W. 

Figure 7 shows good agreement between the data and the Standard Model over the range 
of cos 9w- where there are sufficient data. An important aspect of these results is the 
interrelation between the three spin density correlations ptt, Pll and plt determined in 
this analysis. These quantities can be displayed in one common plot due to the fact that 
they satisfy the condition (7). Thus, they can be plotted in a triangle plot, as shown 
in figure 9. Data can be plotted directly as indicated with arrows on the left side of 
figure 9(a) or by using Cartesian coordinates x and y, 

X = -^(Ptt - Pll) , y ^ Plt ■ 

In figure 9(a) there are seven data points, each one corresponding to a different bin 
of cos 6w-- Error bars are also shown. The curve inside the triangle is the locus of 
points {x,y) calculated using CC03 diagrams. Each point on the curve corresponds to 
one particular value of cos9\y-. The point corresponding to cos9w- = —1 is at the top 



18 




Figure 9: (a) is the joint plot of the TT, LT and LL correlation data shown in figure 7. 
The coordinates x and y are defined in the text. The points are numbered 1-7, the 
first one corresponding to the data bin in cos 6w- at —0.75 and the last one to the bin 
at cos^i^- = +0.9. The curve inside the triangle plot is explained in the text, 
(b) shows the {x,y) point obtained as the average of the points in (a). The shaded area 
is one standard deviation around the average (x, point. The star symbol shows the 
corresponding point predicted by the Standard Model. 

vertex of the triangle. Further points are distributed as indicated by the arrows along 
the curved line and the last point (for cosOw- = +1) is back at the top vertex. The solid 
part of the curve indicates the region where most of our data are located. 

Figure 9(b) shows the average values, ptt, Plt and pll, of the spin density matrix 
elements of Table 1, presented in a triangle plot. The shaded area is the one standard 
deviation region around the average point {x,y). The star symbol is at the point predicted 
by the Standard Model using the CC03 diagrams. 

It is clear that the TT correlation probability is large, the LT correlations are next in 
strength and the LL correlations are small. These are some of the most important features 
of the Standard Model predictions for the reaction e~e~^—>-W~W~^. 

In spite of the limitations due to low statistics, these results show that the salient features 
of the Standard Model predictions for the W — W polarisation correlations are compatible 
with our data. This provides an additional test of the gauge theory relations between the 
SU{2)l ® SU{1)y couplings. 
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